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Abstract 

Biquaternionic Vekua-type equations arising from the factorization of 
linear second order elliptic operators are studied. Some concepts from 
classical pseudoanalytic function theory are generalized onto the consid- 
ered spatial case. The derivative and antiderivative of a spatial pseudoan- 
alytic function are introduced and their applications to the second order 
elliptic equations are considered. 

1 Introduction 

Theory of pseudoanalytic functions [I], [1] is one of the classical branches of 
complex analysis extending the concepts and ideas from analytic function theory 
onto a much more general situation and involving linear elliptic equations and 
systems with variables coefficients. Its development in forties-fifties of the last 
century was fast and deep and historically represented an important impulse to 
the progress in the general theory of elliptic systems. Nevertheless the important 
obstacles for a further development of pseudoanalytic function theory were its 
limited practical applications together with the fact that many important results 
remained in the level of existence without a possibility to make them really 
applicable for solving problems of mathematical physics. Recent progress in 
pseudoanalytic function theory reported in [3] shows its deep relation to the 
stationary Schrodinger equation and more general linear second order elliptic 
equations and includes new results which allow one to make the basic objects of 
the theory and their applications fully explicit. Among other results it is worth 
mentioning the possibility to obtain complete systems of solutions to second 
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order elliptic equations with variable coefEcients and their use for solving related 
boundary and eigenvalue problems |4]. 

The generalization of the concept of a pseudoanalytic function onto the case 
of three or more dimensions was a subject of considerable efforts (see, e.g., [3], 
[S], [5], d], |H])- The development of a spatial generalization of the theory of 
pseudoanalytic functions seems to be possible if only the study is restricted 
to linear first order systems of relatively special forms otherwise it is difficult 
to expect an efficient generalization of the very first concepts of the theory 
like derivative and antiderivative not to mention formal powers, the Cauchy 
integral formula and other notions. In the present work we study the Vekua 
equation arising from the factorization of the stationary Schrodinger operator. 
This biquaternionic Vekua equation was first introduced in [9], [3] as a direct 
generalization of the so-called main Vekua equation closely related to the two- 
dimensional stationary Schrodinger equation. We show that it is possible in 
a natural way to introduce the concept of a derivative for solutions of the bi- 
quaternionic main Vekua equation and exactly as in the two-dimensional case 
the derivative is a solution of another Vekua equation the operator of which is 
involved in the factorization of the stationary Schrodinger operator. Moreover, 
we also show that this operation of differentiation is invertible and introduce 
the notion of the antiderivative. Among other applications these constructions 
allow us to construct exact solutions to the stationary Schrodinger equation 
from the known ones possessing certain symmetries. 

In this work we consider the elliptic case only, however the generalizations 
of the presented results onto the hyperbolic situation (see [10] and [7]) are also 
possible. 

2 Some definitions and results from two-dimensional 
pseudoanalytic function theory 

In this section we give some preliminary definitions and results from the classical 
(complex) pseudoanalytic function theory. For more details, see [T| l3l ITT]. 

Let n, he a. simply connected domain in M^. We say that a pair of complex 
functions F and G possessing in partial derivatives with respect to the real 
variables x and y is a generating pair if it satisfies the inequality 



From this definition it follows that if W : fl ^ <C, then there exist unique real 
valued functions (f>, ip such that 



Im(FG) > 



in ri. 



W{z)^<j){z)F{z)+iP{z)G{z), Vzerj, 



where (p 



Im(M^G) 
lm(FG) 



and Ip 



Im(T4^i^) 
lm(FG) 
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For a given generating pair {F, G) in fi, we can define an {F, G)-derivative 
W{zo) on : ^ C if the (finite) limit 

W{zo) = fim W^(^) - '^(^o)-P^(^) - ^(^o)g(^) 

z~^zo Z — Zq 

exists. Sometimes instead of the notation W for the {F, G)-derivative of W we 

will use the notation ^"^'f — . 

az 

We denote hy dz = ^{dx + idy) and dz = \{dx — idy), where dx '■= ^ and 
dy := The following expressions are known as characteristic coefficients of 
the pair {F, G) 

FGz-FzG ^ FGz-FzG 



a{F,G) — =7r' °iF,G) 



FG-FG ' ' ' FG~FG 



FGz ~^ F^G ^ FGz — FrG 



Af.G) = =7r' ^{F.G) 



FG-FG ' ' ' FG~FG 

where the subindex z or z means the application of or dz respectively. 



Theorem 1 Let {F, G) be a generating pair in and W £ G^{Q). The {F, G)- 
dipQ)W 

derivative W — ' — of W exists and has the form 
dz 

W = (bzF + ,PzG^Wz~A^p^G)W-B^F,G)W (1) 



if and only if 

Wz=a^F,G)W + \f^g)W. (2) 

Equation ([2]) is called Vekua equation (sometimes Carleman-Vekua equation) 
and can also be written in the following form 

^F + ^-G = 0. 

The Vekua equation ^ represents a generalization of the Cauchy-Riemann 
system and solutions of this equation are called (F, G)-pseudoanalytic functions. 

It is easy to see that functions F and G are {F, G)-pseudoanalytic, and F = 

G = 0. 

Let (F, G) and (Fi, Gi) be two generating pairs in {I. We say that (Fi, Gi) 
is a successor of (F, G) and (F, G) is a predecessor of (Fi, Gi) if 

o-(Fi,Gi) = ".{F.G) and ^(_Fi,Gi) = -B{^p^Q). 
The importance of this definition becomes obvious from the following result. 

Theorem 2 Let W he an {F,G)-pseudoanalytic function and let (Fi,Gi) be a 
successor of (F, G). Then W is an {Fi,Gi)-pseudoanalytic function. 
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In other words, when W is an {F, G)-pseudoanalytic function then W satis- 
fies the Vekua equation {W)- = a(^F,G)^ — B^pQ-fW. 

For a given generating pair (F, G) in f2, the adjoint generating pair (F^ G)* = 
{F* ,G*) is defined by the formulas 

F* = and G* 



FG -FG FG- FG 

The {F, G)-integral is defined as follows 

J Wd(^F,G)Z = F{zi)^e J G*Wdz + G{zi)Re J F*Wdz 

where F is a rectifiable curve leading from zo to zi. 

If VF = + V'G is an {F, G')-pseudoanalytic function where (j) and ip are 
real valued functions then 

Wd^F^G)^ = W{z) - </.(zo)F(z) - V(^o)G(z), 



and as F =G = 0, this integral is path-independent and represents the {F, G)- 

antiderivative of W. 

Consider the equation 

^2 = $ 

for a real-valued function 1^9 in a simply connected domain, where $ = $1 + i$2 
is a given complex valued function such that its real part $1 and imaginary part 
^>2 satisfy the equation 

ay$i+9,$2 = 0, (3) 

then we can reconstruct ip up to an arbitrary real constant c in the following 
way 

(p{x, y) = 2 ^idx - <^>2dy^ + c. (4) 
By A we denote the integral operator in (|4]): 



A[<^>]{x,y) = 2(^J <i>idx - <S>2dy 



Thus if $ satisfies ([3]) , there exists a family of real valued functions (p such 
that (fiz = given by the formula (p — A[^]. 

In a similar way we define the operator A corresponding to equation (/s^ = $, 
where $ = $1 -|- satisfies dy^i — dx^2 = 0, by the following antiderivative 
operator: 



^[i'](a;, y) = 2 $1^2: -I- $2^2;^ + c. 
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Now, it is well known that if /o is a non vanishing particular solution of 
the one-dimensional stationary Schrodinger equation ( — + q{x)) f{x) = 0, 

the Schrodinger operator can be factorized as ^ — q = + ^)(^ — 

A similar factorization can be obtained in the two-dimensional case. First, we 

define the Vekua operators 

V-.^ch-^-jC and Vi-.^d^+^jC, 

where / G C'^(ri; M) is a nonvanishing real- valued function defined in Q, and C is 
the complex conjugate operator. The two-dimensional stationary Schrodinger 
equation 

{ - l\ + q{x,y))^^Q inVt, (5) 

where A is the two-dimensional Laplacian A :— -1- ^5-, q and (p are real- 
valued functions, Lp £ C'^{fl;R), can be factorized as follows. 

Theorem 3 [12] Let f be a real-valued nonvanishing solution in of the Schrodinger 
equation Then for any real-valued function ip G C^(57) the following equal- 
ities hold 

^{/\- q)if = ViVip = ViVif. (6) 

Theorem 4 [9] Let W = Wi -\- iW2 he a solution of the Vekua equation VW = 
0. Then Wi — Re W is a solution of the Schrodinger equation in CI and 
W2 = Im W is a solution of the associated Schrodinger equation 

(-A + r)V = mVt, r = 2 (^^^ - g, (7) 

where{Vff^fl + fl. 

Theorem [5] tells us that as much as real and imaginary parts of a complex 
analytic function are harmonic functions, the real and imaginary parts of a so- 
lution of the Vekua equation VW = are solutions of associated stationary 
Schrodinger equations. Now, we know that given an arbitrary real-valued har- 
monic function in a simply connected domain, a conjugate harmonic function 
can be constructed explicitly such that the obtained pair of harmonic functions 
represent the real and imaginary parts of a complex analytic function. This 
corresponds to the more general situation for solutions of associated stationary 
Schrodinger equations; the following theorem gives us the precise result. 

Theorem 5 |12| Let Wi be a real-valued solution of in a simply connected 
domain f7. Then the real-valued function W2, a solution of such that W = 
Wi -\- iW2 is a solution of VW — 0, is constructed according to the formula 

W2^ f-'A[if^d^{f-'W,)]. 

Given a solution W2 of the corresponding solution Wi of (0) such that 
W = Wi -\- iW2 is a solution of VW = 0, is constructed as 

Wi = -fA[tf-^d^{fW2)]. 
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3 Pseudoanalytic derivative for a class of Vekua 
equations in the plane 

For / e C^(f2,M), a nonvanishing real-valued function defined in fl, it is known 
[3J Sect. 3.4] that the functions F — f and G ~ i/ f represent a generating pair 
for the Vekua equation 

VW = 0, (8) 
and if one fixes this generating pair {F, G) = {f,i/f), the equation 

Viw = (9) 

is satisfied by the {F, G')-derivatives of solutions of (|S]) and in particular by the 
successor (Fi, Gi), and hence it is natural to call Vi a successor of V. 
Note that from (HI) the {F, G')-derivative in this case is given as 

W = ^ =W,- = CVGW. 

dz f 

Thus, the fact that is a solution of can be written in the following form. 
If 14^ is a solution of ([5]) then 

ViGVCW = 0. 
Note the following commutation relation 

GV^VG 

where V ^ d,^ - ^G. That is, 

W ^VW 

and ViW — V\VW . The equality to zero of this last expression in the case 
when is a solution of ([8]) is a direct consequence of the factorization ([6]), 
i.e. ViVcp — j{A ~ q)ip where q = ^ and ip is real- valued, twice continuously 
differentiable function, and 

ViVii(j)) = ^(A-r)<^ 

where r = —q + 2 (^^^ : as well as of the fact that from theorem S] if is a 

solution of (HI) then (-A + q)Wi = and (-A + r)W2 = where Wi = ReW 
and W2 = Im W. 

Let us notice that the {F, G')-derivative of W can be also represented as 
follows 

VW = fd.if-'w,) + if-^d^uw^). (10) 
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This form suggests the definition of the inverse operation, the antiderivative: 

w^fA{f-^<i>,) + tf-^A{f^2) 

where <& = <i>i+«$2 is any solution of (|9]). The resulting function w is necessarily 
a solution of ([S]). 

From (|10l) it can be seen that the {F, G)-derivative can be written as 

W = ^,F + i:,G (11) 
where (p — f~^Wi and = /W2. 

4 Three-dimensional pseudoanalytic function 
theory for a class of Vekua equations 

4.1 General definitions and results 

We denote by M{C) ~ {Q \ Q — J2a=oQa^a} the algebra of complex quater- 
nions (biquaternions) (see, e.g., [13l HI]), where Qa G C, eo = 1, {e^ | k = 
1, 2, 3} are the standard quaternionic imaginary units satisfying e^e^ + e^e^ = 
—25a^f3, where Sa,p is the usual Kronecker delta. The imaginary unit in C is 
denoted by i as usual and commutes with all basic units Bq,, a = 0,1,2,3. The 
vector representation of Q S 11(C) can also be used: Q = Sc{Q) + Vec{Q) , where 
Sc{Q) — Qo and Vec{Q) = Q = Qk^k- The quaternionic conjugation is 

defined by Q = Qo — Q (sometimes the quaternionic conjugation operator Ch 
is used, i.e. ChQ — Q)- 

By we denote the operator of multiplication by a biquaternion P from 
the right-hand side, i.e. 

Let Q be a complex quaternion- valued differentiable function of x = (xi, 0:2, X3). 
We define the Dirac operator (sometimes called the Moisil-Theodorescu opera- 
tor) as 

3 

DQ = Y,(^k dkQ, (12) 
where dk ■= t: — . Expression (1121) can be rewritten in a vector form as 

OXk 

DQ = -div Q + VQo + rot Q, 

in other words Sc{DQ) = —div Q and Vec{DQ) = VQo + rot Q. Let us notice 
that = — A, where A := ^^^-^ 9^ is the Laplacian in M.^. 

We notice that the Dirac operator was introduced as acting from the left- 
hand side. The corresponding operator acting from the right-hand side we will 
denote by DrQ — J2k=i ^kQ ^k- In the vector form can be represented as 
DrQ = -div Q + VQo - rot Q. 



7 



The Dirac operator admits the fohowing generaUzation of the Leibniz rule 
(see, e.g., [15]): 

Theorem 6 Let {P,Q} C ]HI(C)) , where D, is some domain inR^. Then 

D[P ■ Q] = D[P] ■ Q + P ■ D[Q] + 2{Sc{PD)) [Q], 

where 

3 

{Sc{PD))[Q]^-Y,PkdkQ- 

k=l 

In particular, we observe that when Vec(P) = 0, i.e. P = Pq, then 

D[Po ■ Q] = D[Po] ■ Q + Po ■ D[Q]. 
Consider now the equation 

where '3/ = '^k^k is a given complex-valued vector such that rot = 0. 

The complex- valued scalar function (p is then said to be the potential of '4'. We 
will write ip — -4['i']. The operator ^ is a simple generalization of the opeators 
A and A of the preceding section. It is well known that we can reconstruct (f, 
up to an arbitrary complex constant c, in the following way 

A[^]{x,y,z) ■i>i{tyo,zo)d^+ ■$2{x,t], zo)dr] + ^'3(0;, y, C)dC + c, 

Jxo Jya J 

where (xq, yo, zq) is an arbitrary point in the domain of interest. 

Let us now consider the three-dimensional stationary Schrodinger equation 

(-A + g)0 = O inn, (13) 

where q and are complex- valued functions depending on x = (a;i,X2,a;3) and 
J7 is a domain in M'^. Using the Dirac operator, one can factorize this equation 
in a way similar to that in theorem [6] 

Theorem 7 [31 Sect. 16.1] Let f be a nonvanishing particular solution of US]) . 
Then for any scalar (complex-valued) function (j) £ C^(Vl] C) the following equal- 
ity holds: 

{-A + q)cf,= {D + M^){D-^CH)<l>. (14) 

The equivalent of theorem |4] in the three-dimensional case is given here. 
Theorem 8 Let Wq + W be a solution of the equation 

D-^Ch]w = 0. (15) 



/ 

Then Wq is a solution of ilS\) with q — the scalar function ipQ — Wq/ f is 
a solution of the equation 

div(/2v^o) - (16) 
and the vector function $ = /W is a solution of the equation 

TOt{f-^ rot *) = 0. (17) 
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4.2 Pseudoanalytic derivative in the space for a class of 
Vekua equations 

We denote the operators 

V:=D-^Ch and Vi:^ D + M^. 

Here the notations from Section |3] are preserved in order to make the analogy 
between the two-dimensional and the higher dimensional case more transparent. 
Application of the operator Ch to gives us 

Df Df 

We denote 

— Df Df 

V:^Dr-M—CH and Vi:^Dr + -^. 

All the operators V, V, V\ and V\ are considered in application to ]HI(C)-valued, 
continuously difFerentiable functions. Note that the following relations hold 

ChV = ~VCh and CrVi = -ViCh- 

Suppose 

VW = Q (18) 
where W G Ci(rj;H(C)). Then due to (HU we have 

FiFM^o = FiFVKo = (-A + g)VKo = (19) 

and hence 

= Fi^W = ChV iVW = ViVW. 

Thus, from this and from (IT^ the H(C)-valued function VW is a solution of 
the equation 

Viw = {) (20) 

whenever W is a, solution of ((T8)) . Consequently we define the Bers derivative 
of a solution W of (fT8|) as follows 

Definition 9 Let W be a solution of hl8\) . Then the Junction 

W = VW 

is a solution of i2U\) and is called the derivative ofW. 

In other words, V : kcrV — ^ ker Vi. Moreover, as we show below the deriva- 
tive W is purely vectorial, and hence it is a solution of the equation Viw — 
as well. 
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Let us obtain another representation for W, similar to (ITT|) . We have 

3 

where represent the generating quartet [31 Sect. 16.2] for equation ((T5)) 
p-f p-^i p_e3 

</?/c are scalar functions defined by (/?o = Wq/ f and (pfc = /M^fe for k = 1,2,3. 
Then 

3 3 3 7— \ p 

Q=0 Q=0 a=0 ■' 



Q=0 a=0 

It is easy to see that VFa = for a = 0, 1, 2, 3. Indeed, we have 

= VF^= ChVF^ = VFa, 

and if a = 0, then VFq ^ VFq = 0, otherwise VFk = -VFk = for fc = 1, 2, 3. 
Thus, another representation for the derivative W has the form 

3 

= = (21) 

a=0 

Let us obtain the inverse operation,- the antiderivative. We have the following 
two equalities 

3 

Y,{D^a)Fn=0 (22) 

and 

3 

Y,iD^c.)F^ = -W, (23) 

0=0 

the first of them being another form of the Vekua equation (|18l) (see [31 Sect. 
16.2]), and the second is obtained from (I?!]) by applying Ch. More explicitly. 



fc=l ■' 

and 

3 



k=l ^ 
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Then _ 

Dfo = -Yf^- (24) 

A scalar solution (^o ol this equation exists if only W is purely vectorial and 
satisfies the equation 

rot {jW^ = 0. 

Consider 

— 3 3 

fe=l •' k=l ^ 

1 ^ 

^ k=l 

= E ^'"'^'^ div * 

fc=i ^ 

where * := J2k=i "fk^k- 

As * = /W we have that [31 equation (16.11)] 



div * = 0. (25) 



Thus, W + W = Q and hence W is purely vectorial. 
Consider the expression 

rot Qiy^ = rot F„ • D^, 



rot ^V(^o + E 



From the previous chain of equalities we already know that 

3 3 

^ GkDLpk = - ^ Dipkek = = - rot *, 

k=l k=l 

because of Thus, 

rot ( \w] = — rot I -^rot4> 



2 
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and according to theorem [5] this is zero (see [31 theorem 161]). 

We have then that if is a solution of ([T8| then equation p4)) possesses a 
solution tpo which can be written as follows 



f 



W 



(26) 



This expression is unique up to an additive constant. 

Now let us see, how $ can be recovered from p2|) and ((23)) . We have 

Q 



k=l 



that is, 



More explicitly. 



f • 

-W. 

2 



div $ = and rot $ 



/ 



(27) 



Notice that div (fW) = 0. Indeed, according to p4|. 

'Wo 



/ 



and due to theorem H div (/^V (^)) = 0. 

Thus, the problem of recovering $ reduces to the well-studied problem of 
reconstruction of a vector by its divergence and rotor. 

We will need the following notation. For a given vector function Q, wc define 
B[Q](x) by 



B[Q](x) 



1 



Q(y) 



Att |x - y| 

Note that B is a right-inverse for the operator rot rot. 
Then we have (see, e.g., [1^1 Sect. 5.7]) 



* = rot B 



V/i, 



(28) 



where h is any harmonic function defined in Vt. 

These results allow us to formulate the following theorem. 



Theorem 10 Let w G C^{^] 
tion 



he a purely vectorial solution of the equa- 
(29) 



D + M"T- ) w 



where f is a nonvanishing scalar -function. Then a solution of the equation 

DW - = inn (30) 



Df— 
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such that Mv = W is defined as follows 



-irot(B [/w]) + ^ 



(31) 



where h is an arbitrary harmonic function in fl. 
Proof. First, wc notice that 



rot ( y I = 0. 



(32) 



Indeed, from the vector part of ([29|) one has rotw + (w x ^) = which is 

equivalent to ((32|). Thus, the vector y is a gradient and the expression A 
is well defined. 

Observe that the operator in equation pop can be written in the following 
form 

Df 1 

D - -j-Cu = 2 {fDf-\l + Ch) + f-^Df{I - Ch)) 



TJ 



where / is the identity operator. Applying it to (1311) we obtain 



{fDf-\l + Ch) + f-'Df{I - Ch)) fA 



1 



f\ f 



rot (B [/w]) 



fDA 



f-'D rot (B [/w] 



(w — w) = 0. 



Next we prove that w — W. For this we observe that 

W^VW^^ {fD.rf-\l + Ch) + f-'Drf{I - Ch)) W 



-2 I ^^-^ 



f-^D iot{B [fw]) = - (w + w) = w. 



Another important corollary of the preceding calculations is the following 
generalization of theorem [5] 

Theorem 11 Let Wq be a scalar solution of U3\} with q ~ ^ in fl. Then 
the vector function W such that # = /W be a solution of ( [j7| ) and 1125]) and 
W = Wq + W be a solution of il5\) . is constructed according to the formula 



W = -/-I {rot (B [/2 V (f-^Wa)] ) + V/i} 



(33) 



Given a solution €> of ( [_?7| j and \25\) . for~W — j'^ the corresponding solution 
Wq of hli^) such that W = VFq+W be a solution of \15\) . is constructed as follows 



Wq - -fA rot(/W)] 



(34) 
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Proof. Let Wo be a scalar solution of (IT^ with q = -j-. Then from we 
have that W = 2/V (/~^Wo) and due to ([28l), we obtain 
* = -rot (B [/' V(/-iT4^o)]) +V/i, 

and as a consequence ((33)) holds. 

Now let * be a solution of ([T7|) and According to we have W = 

-2f-^D^ = - rot (/W). Now from (gH) we obtain ([Ml). ■ 

5 Examples and applications 

In the present section we analize the possibility to construct solutions for the 
equation (1291) under certain restrictions on the function / and apply them for 
obtaining solutions of the Schrodinger equation ([T^ . 

Proposition 12 Let f be a nonvanishing scalar function defined in £7 and 
there exist a harmonic function p such that 

\7f xVp^O in n. (35) 

Then the purely vectorial function 

F = jDp (36) 

is a solution of i29\) and G = fDp is a solution of the equation 

{p-M^^G = Q. (37) 

Proof. Consider DY = -jjDf Dp = -^Dp^ = ^ similar way 

one can check the validity of ([57|. ■ 

Let us notice that (|35|) is fulfilled when / depends on p and this is of a 
considerable importance because harmonic functions are behind all physically 
interesting orthogonal coordinate systems. For example, consider the spherical 
coordinate r = x'^ + y'^ + and let / = f(r). The function r itself is not 
harmonic, however the function p = 1/r is harmonic everywhere except the 
origin. Thus, when / is a function of r and hence a function of p, and the 
domain f2 does not contain the origin the function (fSG)) is a solution of ([29]) . 

Proposition 13 Let f be a nonvanishing scalar function defined in and 
there exist a harmonic function p such that 

(V/, Vp) =0 mn. 

Then the purely vectorial function 

F = fDp 

is a solution of i29\) and G = jDp is a solution of {5^. 
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Proof. Consider DF = Df Dp = —fDp^ = — F-^. In a similar way one 
can check the vahdity of ([37| . ■ 

These two propositions give us the possibihty to obtain three independent 
solutions of equation ((29)) in a simple way when / is a function of a harmonic 
function and this harmonic function possesses two different orthogonal harmonic 
functions. As we show in the next example this situation is of practical impor- 
tance. Let us notice that in general for a given harmonic function in three dimen- 
sions there not necessarilly exists a nonconstant orthogonal harmonic function 
(a counterexample is given in jl7|). 



Example 14 Consider the cylindrical coordinate system r — \J + y^, = 
arctan(2//a;) and z. Suppose that f depends only on r. Then the required har- 
monic function p can be chosen as p = logr, and obviously the functions 9 and 
z are harmonic and orthogonal to p. Thus Fi = j^Dp = (^^i ~^ '^^2) , 

F2 = f{r)D9 = -f{r) (-^ei + ^ea) and F3 ^ /(r)e3 are solutions of ^ 
independent in the sense that any function w defined in a domain not contain- 
ing the axis r = and having purely quaternionic values can he represented in 
the form w = Vk^k where ip^ are scalar functions. Consequently Fi, F2 

and F3 form a generating triplet of equation i f^ff)) ( see more on this concept in 
^) and w is a solution of i29\) if and only if the scalar functions (pk satisfy the 
equation 

3 

Y^DifikFk^O mil. (38) 

fe=i 

In general we note that if f is a function of a harmonic function depending 
on two Cartesian variables then we always can obtain a generating triplet for 
equation h29\l and hence reduce the equation to the form ifg^|J. Indeed, for such 
harmonic function we always have two orthogonal harmonic functions: its con- 
jugate harmonic depending on two Cartesian variables as well, and the third 
Cartesian coordinate. 

Assume that / is a function of a harmonic function p. Then according to 
proposition [121 F ~ J^P a- solution of ([29]). Its antiderivative ([31]) gives us 
a solution of (pO| and according to theorem [8] the scalar part of the obtained 



solution is a solution of (|T3)) . Thus we obtain that -0 = f-^ 



IS 



necessarily a solution of (|T5| . and hence the results of this section allow us to 
construct new solutions of the Schrodinger equation (I13p from known solutions 
possessing certain symmetries. 



6 Conclusions 

For biquaternionic Vekua type equations arising from the factorization of linear 
second order elliptic operators we introduced and studied the concepts of the 
derivative and of the antiderivative in the sense of Bers. These concepts rep- 
resent new relations between first order elliptic systems and allow us to study 
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solutions of the second order elliptic equations in relation to components of 
the biquaternion valued solutions of biquaternionic Vekua equations. We have 
shown that the revealed relations are useful for obtaining new exact solutions of 
the Schrodinger type equations and at the same time they imply natural open 
questions important for the general theory of linear second order elliptic equa- 
tions. Namely, we showed how in a quite general and important for practical 
applications situation one can obtain a generating triplet for the biquaternionic 
Vekua equation which can be considered as a successor of the main Vekua equa- 
tion (15(11) . In the two-dimensional case this was sufficient to obtain a whole 
generating sequence related to the Vekua equation and as a consequence a com- 
plete system of solutions called formal powers (see [3]) for the Vekua equation 
and for the related second order elliptic equation. At present it is not clear how 
to generalize those results onto the spatial situation. 
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